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ABSTRACT 

While O is often seen in spectra of Type la supernovae (SNe la) as both unburned 
fuel and a product of C burning, C is only occasionally seen at the earliest times, and 
it represents the most direct way of investigating primordial white dwarf material and 
its relation to SN la explosion scenarios and mechanisms. In this paper, we search 
for C absorption features in 188 optical spectra of 144 low-redshift (z < 0.1) SNe la 
with ages ^3.6 d after maximum brightness. These data were obtained as part of the 
Berkeley SN la Program (BSNIP; Silverman et al. 2012a, b, submitted) and represent 
the largest set of SNe la in which C has ever been searched. We find that ~11 per cent 
of the SNe studied show definite C absorption features while ^^25 per cent show some 
evidence for C II in their spectra. Also, if one obtains a spectrum at t < —5 d, then 
there is a better than 30 per cent chance of detecting a distinct absorption feature 
from C II. SNe la that show C are found to resemble those without C in many respects, 
but objects with C tend to have bluer optical colours than those without. The typical 
expansion velocity of the C II A6580 feature is measured to be 12,000-13,000 km s^^, 
and the ratio of the C II A6580 to Si II A6355 velocities is remarkably constant with 
time and among different objects with a median value of ~1.05. While the pseudo- 
equivalent widths (pEWs) of the C II A6580 and C II A7234 features are found mostly 
to decrease with time, we see evidence of a significant increase in pEW between ~12 
and lid before maximum brightness, which is actually predicted by some theoretical 
models. The range of pEWs measured from the BSNIP data implies a range of C mass 
in SN la ejecta of about (2-30) xlO"^ M©. 
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1 INTRODUCTION 



It is thought that thermonuclear explosions of C/0 white 
dwarfs (WDs) give rise to Type la supernovae (SNe la; e.g., 
Hoyle fc Fowler|1960| [Colgate fc McKce.1969, Nomoto et al.| 



1984[ see [Hillebrandt fc Niemeyer||2000| for a review). How- 
ever, after decades of observations and theoretical work, the 
details of SN la progenitors and explosion mechanisms are 
still missing. Despite this, SNe la have been used in the re- 
cent past to discover the accelerating expansion of the Uni- 
verse ( Riess et al.||1998 Perlmutter et al.|[1999 l, as well as 



to measure cosmological parameters (e.g., 
Riess et al. 2007 



Astier et al. 



licken et al 



2006 



2009 



Wood-Vasey et al.|2007[ 

Kessler et "aLl|2d09|7Amanullah et al.||2010[ [Suzuki eTlil 
20T2I 



As the explosion proceeds, pristine C and O from the 
progenitor WD may get mixed throughout various layers of 
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the ejecta. However, the amount and exact location of this 
unburned material varies widely among published models 
(e.g., ,H6flich et al.||2002| [Gamezo et al.||2003| |R6pke et al 



2007 



Kasen et al.|2009 l. Thus, determinations of the quan- 



tity of these elements after explosion, their spatial distribu- 
tion in the ejecta, and other measurements of this primordial 
material can help constrain possible explosion mechanisms 
of SNe la. 

Oxygen is found in SNe la as both unburned fuel from 
the progenitor WD as well as burned ash as a product of 
C burning. Therefore, it is difficult to associate O detec- 
tions, which are common in SNe la (e.g., |Filippenko|[T997l ), 
with primordial material. This leaves us with C as the most 
direct link to matter from the pre-explosion WD. Observa- 
tions during the first 2-3 weeks after explosion probe the 
outermost layers of the ejecta, which is where the unburned 
material is most likely to reside. At the temperatures and 
densities observed in SNe la at these epochs, singly ionised 
states of C are expected to be the dominant species (e.g.. 
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Tanaka et al.|2008 L Neutral C could appear (mostly in the 
near-infrared) at significantly lower temperatures, and dou- 



bly ionised C would require much higher temperatures ( Mar- 



ion et al. 20061. Thus, the best chance of detecting unburned 



material is to look for C II features before and near _B-band 
maximum brightness. 

Until recently, C detections in SNe la were only no- 
ticed in a small handful of objects. Most of the SNe la that 
show obvious C II absorption features are extremely lumi- 
nous objects with slowly evolving light curves and excep- 
tionally low expansion velocities. These objects are thought 



to arise from super-Chandrasekhar-mass WDs ( Howell et al. 
2006i [Yainanaka et al.||2009| |Scalzo et al.||2010| |Silverman 



et al. 2011 



Taubenberger et al. 20111, and thus the detec- 



tion of unburned C in their spectra is likely related to the 
relatively rare explosion mechanism that produces these ob- 
jects. In addition, there are a few instances of relatively nor- 
mal SNe la (i.e., ones that obey the light-curve decline rate 
versus luminosity at peak brightness relation, known as the 
"Phillips relation" ; Phillips|1993 1 that also show strong C II 
absorption (e.g., Patat et al.|1996 Garavini et al.|2005 l. On 
the other hand, it has usually been found that C does not 
appear at all in early-time optical spectra of SNe la, or it is 
weak or extremely blended (e.g., Mazzali|2001| Branch et al 



20031 IStanishev et al.||2007| iThomas et al. 20071 



Recently, however, C detections in SNc la at early times 
have become more common thanks to the amassing of many 
more spectra of SNe la at very early epochs, higher signal-to- 
noise ratio (S/N) data, and astronomers being more meticu- 
lous in their search for this elusive unburned material. IPar-l 
[rent et ah] ( |2011[ ) present new spectra of 3 SNe la that show 
distinct C absorptions and analyze those alongside 65 other 
objects from the literature. They find that C II features are 
detected more often than previously thought and estimate 
that 30 per cent of all SNe la may show evidence of unburned 
C in their pre-maximum spectra. This work is explained by 



Thomas et al. (20111, who discuss observations and analy- 



ses of 5 more objects that show C II features. They conclude 
that the C is likely distributed spherically symmetrically and 
that per cent of SNe la show C absorption features 

at epochs near 5 d before maximum brightness. Finally, |Fo-| 



latelli et al. (20121 use data from the Carnegie Supernova 



Project to determine that at least 30 per cent of objects 
show C II absorption and that the mass of C in the ejecta 
is consistent with 10~^-10~^ Mq. They also find evidence 
that SNe la with C tend to have bluer colours and lower 
luminosities at maximum light. 

In this work, we search for possible C signatures in low- 
redshift (2 < 0.1) optical spectra of SNe la obtained as 
part of the Berkeley SN l a Program (BSNIP). T he data 
are presented in BSNIP I |Silverman et al. 20121, and we 
utilise the spectral feature measurement tools described in 
BSNIP II ( [Silverman, Kong fc Filippenko|2012t . With such 
a large, self-consistent dataset, we are able to accurately 
explore the incidence rate of C in SN la spectra and how 
it varies as a function of observed epoch. Furthermore, we 
can quantify the amount and location in the ejecta of the 
unburned C by measuring pseudo-equivalent widths (pEWs) 
and expansion velocities, respectively. 

The spectral and photometric data used herein are sum- 
marised in Section [2] and our methods for determining the 
presence or absence of C and (for SNe la with definite C 



detection) measuring C II spectral features is described in 
Section|3] Section |4]presents the rate of C detection, the dif- 
ferences and similarities between SNe la with and without C, 
and a discussion of the C II spectral feature measurements. 
We present our conclusions in Section |5] 



2 DATASET 

The SN la spectral data investigated in the current study are 
a subset of those used in BSNIP II and originally published 
in BSNIP I. The majority of the spectra were obtained using 
the Shane 3 m telescope at Lick Observatory with the Kast 
double spectrograph ( [Miller fc Stone|1993| ), and the typical 
wavelength coverage is 3300-10,400 A with resolutions of ~ 
11 and ~6 A on the red and blue sides (crossover wavelength 
~5500 A), respectively. For more information regarding the 
observations and data reduction, see BSNIP I. 

In BSNIP II, we ignored a priori the extremely peculiar 
SN 2000CX (e.g., [Li et al.[[2001[ ), SN 2002 cx (e.g.,|Li et al 



2003 



2006 



2009 



Jha et al. 



20061, SN 2005hk (e.g., Chornock et al 



Phillips et al.|20b7), and SN 2008ha (e.g., [Foley et al 



Valenti et al.||2009j . This was mainly due to the fact 



that they are so spectroscopically distinct from the bulk of 
the SN la population that their spectral features are diffi- 
cult to measure in the same way as for the other objects. 
It should be noted that IParrent et al. (20111 find that all 



of these objects show evidence for unburned C. However, in 
this work we will only concentrate on SNe la that follow the 
Phillips relation, and thus can be used as cosmological dis- 
tance indicators. This means that we also remove all super- 
Chandrasekhar-mass SNe la from our sample, even though 
they show strong absorption from C II (as mentioned above). 

BSNIP II contains 432 spectra of 261 SNe la with ages 
younger than 20 d (rest frame) past maximum brightness. 
We began the current study by inspecting all 206 spectra (of 
156 objects) younger than 5 d past maximum for possible 
C signatures. The oldest spectrum to show evidence of a 
C II feature was obtained ~3.6 d after maximum. Hence, 
for the rest of this study, we only consider spectra younger 
than this epoch. This yields a sample of 188 spectra of 144 
SNe la, which is the largest set of SNe la that has ever been 
inspected for C features. A summary of these objects, their 
"Carbon Classification" (see Section [3| , and their spectral 
classifications based on various classification schemes can 
be found in Table [l] For comparison, [Parrent et aL (20111 
investigate d 58 object j'] y ounger than 1 d past maximum 
brightness, Thomas et al. (2011 1 used 124 objects at epochs 



before 2.5 d past maximum, and the study by |Folatelli et aL\ 
(20121 utilised 51 SNe la with spectra before maximum. 

The spectral ages of the BSNIP data referred to 
throughout this work are calculated using the redshift and 
Julian Date of B-band maximum brightness presented in 
Table 1 of BSNIP I. Furthermore, photometric parameters 
(such as light-curve width and colour information) used in 
the present study can be found in Ganeshalingam et al. (in 
preparation). 



^ When only including SNe la that follow the Phillips relation. 
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Table 1: Summary of SNe la in the Sample 



Namp 


C TT 


SNID 


JJcllt; h hi 


iJi dllC-li 


V V "liJ^^ 


SN Name 


C TT 

\^ 11 


Oi> lU 


Xjcllc L LI 


XjI cLllCll 


V V dligi 




Type 


(bub) iype 


Type 


Type 


Type 






„ ^a 

Type 


(bub)Type 


Type 


Type 


Type 


bJN 1994D 


A 


la-norm 


LVG 


CN 


N 


SN 1998dm 


A 

A 


la-norm 








bN iOOicr 


A 


la-norm 








SN 2002hw 


A 

A 


la-norm 








OAT or\noi c 

SN 2003kr 


A 

A 


la-norm 








SN 2004ey 


A 

A 


la-norm 








OAT or\nr r 

bN iOOocr 


A 

A 


la-norm 


HVG 


CN 


N 


SN 2005el 


A 

A 


la-norm 


LVG 


CN 


N 


OAT tr 

bN iOOoeu 


A 

A 


la-norm 








SN 2005iq 


A 

A 


la-norm 








OAT i^/A/^r'l •] 

SN 2005ki 


A 


la-norm 


LVG 


BL 


N 


SN 2007F 


A 


la-norm 




SS 


N 


OAT (tr\f\'-7 

bN 2007aP 


A 


la-norm 


HVG 


BL 


N 


SN 2007bm 


A 


la-norm 








SN 2007cq 


A 


la 








SN 2008sl'' 


A 


la-norm 




BL 


N 


OAT 1 r\r\ r tt> 

bN 199bbj 


F 


la-norm 




CN 


N 


SN 


1998dk 


F 


la-norm 






HV 


SN 1999dk 


F 


la-norm 








SN 


2000dn 


F 


la-norm 


LVG 


BL 


N 


SN 200 laz 


F 


la-norm 






N 


SN 


2001cp 


F 


la-norm 






N 


OAT i~ir\i^ 1 r 

SN 2001fe 


F 


la-norm 




SS 


N 


SN 


2002ck 


F 


la-norm 




CN 


N 


SN 2002er 


F 


la-norm 






HV 


SN 


2003U 


F 


la-norm 




BL 


N 


SN 2004fz 


F 


la-norm 








SN 


2005ao 


F 


la 




SS 




OAT (tr\f\f 

SN 2005na 


F 


la-norm 


HVG 




N 


SN 


2006ax 


F 


la-norm 








OAT ^^/^f^/^l i 

bN 2006bt 


F 


la-norm 


HVG 


CL 


N 


SN 


2006CS 


F 


Ia-91bg 




CL 




SN 2007bc 


F 


la-norm 




CL 


N 


SN 


2007on 


F 


la-norm 




CL 


N 


SN 2008Z 


F 


Ia-99aa 








SN 


2008hs 


F 


la-norm 








SN 1989M 


N 


la-norm 


HVG 


BL 


HV 


SN 


1991bg 


N 


Ia-91bg 


FAINT 






SN 1994S 


N 


la-norm 




SS 


N 


SN 


1997Y 


N 


la-norm 




BL 


N 


SN 1997br' 


N 


Ia-91T 








SN 


1997do 


N 


la-norm 








SN 1998ef 


N 


la-norm 








SN 


1998es 


N 


Ia-99aa 




SS 




SN 1999aa 


N 


Ia-99aa 




SS 




SN 


1999ac 


N 


la-norm 


HVG 


CN 


N 


SN 1999da 


N 


Ia-91bg 


FAINT 


CL 




SN 


1999dq 


N 


Ia-99aa 


HVG 


SS 




SN 1999gd 


N 


la-norm 




BL 


N 


SN 


2000cp 


N 


la-norm 






N 


SN 2000dg 


N 


la-norm 




SS 


N 


SN 


2000dk 


N 


la-norm 


FAINT 


CL 


N 


SN 2000dm 


N 


la-norm 


HVG 


BL 


N 


SN 


2001bf 


N 


la-norm 






N 


SN 2001br 


N 


la-norm 




BL 


HV 


SN 


2001da 


N 


la-norm 


HVG 


BL 


N 


SN 200 leh 


N 


Ia-99aa 




SS 




SN 


2001ep 


N 


la-norm 


HVG 


CL 


N 


SN 200 lex 


N 


Ia-91bg 








SN 


2002aw 


N 


la-norm 






N 


SN 2002bf 


N 


la-norm 




BL 


HV 


SN 


2002bo 


N 


la-norm 


HVG 




HV 


SN 2002cf 


N 


Ia-91bg 




CL 




SN 


2002CS 


N 


la-norm 








SN 2002cu 


N 


la-norm 








SN 


2002dj 


N 


la-norm 








SN 2002dk 


N 


Ia-91bg 




CL 




SN 


2002eb 


N 


la-norm 






N 


SN 2002eu 


N 


la-norm 


HVG? 


CL 


N 


SN 


2002fb 


N 


Ia-91bg 




CL 




SN 20021ia 


N 


la-norm 


LVG 


BL 


N 


SN 


2002he 


N 


la-norm 


HVG 


BL 


HV 


SN 2002hu 


N 


Ia-99aa 








SN 


2003W 


N 


la-norm 








SN 2003cq 


N 


la-norm 






HV 


SN 


2003gt 


N 


la-norm 








SN 20031ie 


N 


la-norm 


LVG 


BL 


N 


SN 


2003iv 


N 


la-norm 


HVG? 


CL 


N 


SN 2004as 


N 


la-norm 




BL 


HV 


SN 


2004br 


N 


la-norm 








SN 2004bv 


N 


Ia-91T 








SN 


2004bw 


N 


la-norm 






N 


SN 2004dt 


N 


la-norm 


HVG 


BL 


HV 


SN 


2004ef 


N 


la-norm 






HV 


SN 2004eo 


N 


la-norm 








SN 


2004fu 


N 


la-norm 


HVG? 


BL 


HV 


SN 2004gs 


N 


la-norm 




CL 


N 


SN 


2005W 


N 


la-norm 




BL 


N 


SN 2005agj 


N 


la-norm 






N 


SN 


2005bc 


N 


la-norm 


LVG 


CL 


N 


SN 2005de 


N 


la-norm 


HVG 


BL 


N 


SN 


2005er 


N 


Ia-91bg 


HVG? 


CL 




SN 2005eq 


N 


Ia-99aa 


HVG 






SN 


20051z 


N 


la-norm 




BL 


N 


SN 900'iTnei 


N 


Tsi— nrivm 

±(X1HJ1 111 


HVG 




N 


SN 


2006N 


N 


J-d llUl 111 


HVG 


BL 


N 


SN 2006S 


N 


la-norm 


LVG 






SN 


2006X 


N 


la-norm 




BL 


HV 


SN 2006bz 


N 


Ia-91bg 




CL 




SN 


2006cj 


N 


la-norm 




SS 


N 


SN 2006cp 


N 


la-norm 








SN 


2006cq 


N 


la-norm 






N 


SN 2006ef 


N 


la-norm 




BL 


HV 


SN 


2006ej 


N 


la-norm 


HVG 


BL 


HV 


SN 2006et 


N 


la-norm 


LVG 


SS 


N 


SN 


2006gt 


N 


Ia-91bg 








SN 2006ke 


N 


Ia-91bg 


HVG? 






SN 


2006kf 


N 


la-norm 




CL 


N 


SN 20061e 


N 


la-norm 








SN 


2006or 


N 


la-norm 




BL 


N 


SN 2006sr 


N 


la-norm 


HVG 


BL 


HV 


SN 


2007A 


N 


la-norm 


LVG? 


CN 


N 


SN 2007N 


N 


la 




CL 




SN 


2007O 


N 


la-norm 




SS 


N 
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Table 1 — Continued 
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PT 
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9(10 7fr 
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PT 
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N 
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RT 
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SN 2007qe 


N 


la-norm 






HV 


SN 2008ar 


N 


la-norm 




CN 


N 


SN 2008ec 


N 


la-norm 


LVG 


CL 


N 


SN 2008ei 


N 


la-norm 


HVG* 


BL 


HV 


SN 2008s5' 


N 


la 


LVG* 


















SN 2000fa 


? 


la-norm 






N 


SN 2002cd 


7 


la-norm 


LVG 




HV 


SN 2002de 


? 


la-norm 


HVG 


CL 


HV 


SN 2003Y 


7 


Ia-91bg 








SN 2003gn 


7 


la-norm 








SN 2005dv 


7 


la-norm 




BL 


HV 


SN 2006cm 


7 


la-norm 






N 


SN 2006CZ 


7 


Ia-99aa 








SN 2006gr 


7 


la-norm 








SN 20061f 


7 


la-norm 








SN 2007al 


7 


Ia-91bg 








SN 2008bt 


7 


Ia-91bg 




CL 




SN 2008dx 


7 


Ia-91bg 


FAINT* 


CL 

















^Classification based on the presence or absence of C II absorption in our spectra. 'A' — clearly separated spectral absorption 
feature attributed to C II is present; 'F' = a flattening or depression in the red wing of the Si II A6355 feature is present 
that is likely due to C II; 'N' = no evidence for C II is present; '?' = no definitive determination can be made about the 
presence or absence of C II owing to noisy data. 

''Spectral classification using the SuperNova IDentification code (SNID; Blondin fc Tonry||2007 l taken from Section 5 of 
BSNIP I. 



'^Classification based on the velocity gradient of the Si II A6355 line (^Benetti e t al. [2005 1. 'HVG' — high velocity gradient; 
'LVG' — low velocity gradient; 'FAINT' = faint/underluminous. Classifications marked with a '?' are uncertain since light- 



curve shape information is unavailable. Classifications marked with a '*' use the MLCS2k2 A parameter ( Jha et al.|2007 1 
as a proxy for Amis. Taken from BSNIP II 



Classification based on the pseudo-equivalent widths of the Si II A6355 and Si II A5972 lines ( Branch et al. |2009t . 'CN' = 
core normal; 'BL' — broad line; 'CL' — cool; 'SS' — shallo w silicon. Taken fro m BSNIP II. 

'^Classification based on the velocity of the Si II A6355 line ( Wang et al.|2009 l. 'HV = high velocity; 'N' = normal. Taken 



from BSNIP II 

'C II classification was changed from 'N' based on the data presented by 
^C II classification was changed from 
''Also known as SNF20080514-002. 
'C II classification was changed from 

^C II classification was changed from '?' based on the data presented by 
■"Also known as SNF20071021-000. 
'Also known as SNF20080909-030. 



'N' based on the data presented by 
?' based on the data presented by 



Thomas et al. 



Farrent et al. 



(20111 



(21JTT 



Parrent et al. 


2on\ 


Folatelli et al. 


(2012 
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3 CARBON DETECTION AND 
MEASUREMENT 

As mentioned above, C II is the dominant species of C 
for typical SN la temperatures before and near maximum 



brightness (~10,000 K; e.g., Hatano et al. 1999 1. This species 



has four major absorption lines in the optical regime: A4267, 
A4745 , A6580, and A7234 (e.g., |Mazzah|pOOl] [Branch et al.| 
2003 1 . The bluest two lines are usually overwhelmed in SN la 



spectra by broad, blended absorption from iron-group ele- 
ments (IGEs), so we do not attempt to search for either of 
them in our data. The A7234 line is more promising, but still 
not perfect since it is relatively weak, usually quite broad, 
and ofte n falls close to the telluric absorption feature at 
6900 A (jFolatelh et al.||2012|. When this line is detected. 



however, it is often only observed as an "inflection" in the 
spectral continuum ( [Thomas et al.|201ip , though it becomes 
more obvious when the A6580 line is easily detected ([Parrent 
et al.|2011| l 



The A6580 line is the most obvious C II absorption fea- 
ture in the optical (e.g., Hatano et al.[[T999 l and represents 
the best chance of making a definitive detection of unburned 
C in pre-maximum optical spectra of SNe la. However, even 
though it is the deepest C II absorption, it is still signifi- 
cantly weaker than the nearby Si II A6355 line. Furthermore, 
C II A6580 is usually blueshifted to ~6300 A, which often 
intersects the red wing (or emission component of a P-Cygni 
profile) of the Si II A6355 line. Therefore, even though we 
concentrate mainly on searching for C II A6580, unambigu- 
ously observing this feature is still a difficult task. 



help identify spectral features seen in SNe. We fit all spec- 
tra initially classified as 'A' or 'F' using SYNOW, both with 
and without C II, to investigate whether the addition of 
C II significantly improves the match between the synthetic 
and observed spectra. Again, this is similar to previous work 
( Parrent et"al]|2011[ [Folatelh et~al][2012[ ), though [Thomas 



et al. (2011 1 use a different spectral synthesis code, SYNAPPS 



( Thomas, Nugent fc Meza 2011) . 

After using SYNOW to fit all spectra thought to have C, 
it was found that the addition of C II did not improve the 
fit to 33 spectra that were initially classified as 'F' (and 
thus they were reclassified as 'N'). This likely implies that 
the initial visual inspections were perhaps a bit too "op- 
timistic" in detecting possible C absorptions. However, all 
spectra initially classified as 'A' were confirmed to contain 
C in their spectra from the SYNOW fits. Examples of observed 
and synthetic spectra of an A' spectrum, an 'F' spectrum, 
and an 'N' spectrum can be found in the top, middle, and 
bottom panels of Figure [l] respectively. Furthermore, Ta- 
ble [2] lists each spectrum in our sample, along with its age 
and carbon classification, and a summary of the number of 
spectra in each class is presented in Table [3] As mentioned 
above, the oldest spectrum that shows evidence for C II ab- 
sorption (i.e., an 'F' classification) is ~3.6 d after maximum 
brightness. The oldest spectrum with an 'A' classification 
was obtained about 4.4 d before maximum. 



3.1 The Search for Carbon 

The first step in our search for C consists of visually inspect- 
ing each of the 188 spectra in our sample. We concentrate on 
the region 5600-7800 A in order to cover the spectral range 
around Si II A6355 and C II A6580 as well as C II A7234 and 
the O I triplet (centered near 7770 A). Spectra where there is 
an obvious, distinct absorption feature likely associated with 
C II A6580 is given an A' ("Absorption") classification. 'A' 
spectra often also show depressions or distinct absorption 
features associated with C II A7234. Spectra that show no 
distinct absorption, but the possibility of a depression or 
flattening of the red wing of the Si II A6355 feature are clas- 
sified as 'F' ("Flattened"). These data represent tentative 
C detections, and no 'F' spectra show obvious evidence for 
C II A7234 absorption. We consider any spectrum with an 
A' or 'F' classification to "have C" or be "C positive." 

'N' ("No C") classifications are given to spectra where 
there is no evidence of C absorption features and the red half 
of the Si II A6355 line appears to be unaffected by any other 
species. Finally, spectra where no definite classification can 
be made (often due to low S/N; i.e., the noise fluctuations 
were as large as possible C absorption features) are denoted 
as '?' ("Inconclusive"). This classiflcation scheme is similar 
to those used previously ( Parrent et al.|2011 Folatelli et al 



[20T2I ). 

To more quantitatively determine a spectrum's C clas- 



siflcation, we use the spectrum-synthesis code SYNOW ( Fisher 



et al. 19971. SYNOW is a parametrised resonance-scattering 
code which allows for the adjustment of chemical composi- 
tion, optical depths, temperatures, and velocities in order to 
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Figure 1. Examples of observed and synthetic spectra (created 
using SYNOW). (Top) SN 2005iq at 5.9 d before maximum bright- 
ness is an 'A' spectrum, and the synthetic spectrum is shown with 
and without C II included (note that the two synthetic spectra 
are nearly identical at wavelengths below ~6200 A). (Middle) 
SN 2008hs at 7.9 d before maximum is an 'F' spectrum, and 
again the synthetic spectrum is shown with and without C II. 
{Bottom) SN 2008ar at 8.9 d before maximum is an 'N' spec- 
trum, and thus no C II is included in the synthetic spectrum. 
The data are all deredshifted and dereddened using the redshift 
and reddening values presented in Table 1 of BSNIP I, assuming 

that the extinction follows the ICardelli et al.l JTflSSl extinction nni n td a O a jmo a a r»r»r» mrnTTl 
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Table 2: Summary of Investigated Spectra 



SN Name 


Phase'' 


C II 

Type"" 


SN Name 


Phase'' 


C II 

Type'' 


SN Name 


Phase" 


C II 

Type'' 


SN 1989M 


2.49 


N 


SN 


1989M 


3.48 


N 


SN 


1991bg 


0.14 


N 


SN 199 Ibg 


1.14 


N 


SN 


1994D 


-12.31 


A 


SN 


1994D 


-11.31 


A 


SN 1994D 


-9.32 


A 


SN 


1994D 


-7.67 


A 


SN 


1994D 


-5.32 


A 


SN 1994D 


-3.87 


F 


SN 


1994D 


-3.33 


N 


SN 


1994S 


1.11 


N 


SN 1995E 


-2.46 


F 


SN 


1997Y 


1.27 


N 


SN 


1997br 


-4.84 


? 


SN 1997do 


-5.67 


N 


SN 


1998dk 


-7.24 


F 


SN 


1998dk 


-0.54 


N 


SN 1998dm 


-12.48 


A 


SN 


1998dm 


-5.61 


F 


SN 


1998ef 


-8.62 


N 


SN 1998es 


0.28 


N 


SN 


1999aa 


0.24 


N 


SN 


1999ac 


-3.70 


N 


SN 1999ac 


-0.89 


N 


SN 


1999da 


-2.12 


N 


SN 


1999dk 


-6.60 


F 


SN 1999dq 


-3.93 


? 


SN 


1999dq 


2.97 


N 


SN 


1999gd 


-1.12 


N 


SN 2000cp 


2.92 


N 


SN 


2000dg 


-5.09 


N 


SN 


2000dk 


1.00 


N 


SN 2000dm 


-1.63 


N 


SN 


2000dn 


-0.94 


F 


SN 


2000fa 


-8.25 


? 


SN 200 laz 


-3.24 


F 


SN 


2001bf 


1.22 


N 


SN 


2001br 


3.47 


? 


SN 2001br 


3.48 


N 


SN 


2001cp 


1.39 


F 


SN 


2001da 


-1.12 


N 


SN 200 leh 


3.26 


N 


SN 


2001ep 


2.83 


N 


SN 


2001ex 


-1.82 


N 


SN 200 Ife 


-0.99 


F 


SN 


2002aw 


2.10 


N 


SN 


2002bf 


2.97 


N 


SN 2002bo 


-11.94 


N 


SN 


2002bo 


-1.08 


N 


SN 


2002cd 


1.10 


? 


SN 2002cf 


-0.75 


N 


SN 


2002ck 


3.64 


F 


SN 


2002cr 


-6.78 


A 


SN 2002cs 


-7.76 


N 


SN 


2002CU 


-5.28 


N 


SN 


2002de 


-0.32 


? 


SN 2002dj 


-7.98 


N 


SN 


2002dk 


-1.23 


N 


SN 


2002eb 


1.68 


N 


SN 2002er 


-4.58 


F 


SN 


2002eu 


-0.06 


N 


SN 


2002fb 


0.98 


N 


SN 2002ha 


-0.85 


N 


SN 


2002he 


-1.03 


N 


SN 


2002he 


0.29 


N 


SN 2002he 


3.22 


N 


SN 


2002hu 


-5.81 


N 


SN 


2002hw 


-6.27 


A 


SN 2003U 


-2.55 


F 


SN 


2003W 


-5.06 


N 


SN 


2003Y 


-1.74 


? 


SN 2003cq 


-0.15 


N 


SN 


2003gn 


-5.38 


? 


SN 


2003gt 


-5.07 


N 


SN 2003he 


2.71 


N 


SN 


2003iv 


1.76 


N 


SN 


2003kf 


-7.50 


A 


SN 2004as 


-4.36 


N 


SN 


2004br 


3.50 


N 


SN 


2004bv 


-7.06 


N 


SN 2004bw 


-10.03 


N 


SN 


2004dt 


-6.46 


? 


SN 


2004dt 


1.38 


N 


SN 2004ef 


-5.52 


N 


SN 


2004eo 


-5.57 


N 


SN 


2004ey 


-7.58 


A 


SN 2004fu 


-2.65 


N 


SN 


2004fu 


2.43 


N 


SN 


2004fz 


-5.18 


F 


SN 2004gs 


0.44 


N 


SN 


2005W 


0.59 


N 


SN 


2005ag 


0.53 


V 


SN 2005ao 


-1.29 


F 


SN 


2005ao 


0.52 


N 


SN 


2005bc 


1.55 


N 


SN 2005cf 


-10.94 


A 


SN 


2005cf 


-2.11 


F 


SN 


2005cf 


-1.19 


N 


SN 2005de 


-0.75 


N 


SN 


2005dv 


-0.57 


? 


SN 


2005el 


-6.70 


A 


SN 2005el 


1.22 


F 


SN 


2005er 


-0.26 


? 


SN 


2005er 


1.67 


N 


SN 2005eq 


-2.98 


N 


SN 


2005eq 


0.66 


? 


SN 


2005eu 


-9.06 


A 


SN 2005eu 


-5.46 


F 


SN 


2005iq 


-5.86 


A 


SN 


2005ki 


1.62 


N 


SN 20051z 


0.58 


N 


SN 


2005ms 


-1.88 


N 


SN 


2005na 


0.03 


F 


SN 2005na 


1.03 


F 


SN 


2006N 


-1.89 


N 


SN 


2006N 


-0.90 


N 


SN 2006S 


-3.93 


? 


SN 


2006S 


2.99 


N 


SN 


2006X 


3.15 


N 


SN 2006ax 


-10.07 


F 


SN 


2006bt 


-5.30 


F 


SN 


2006bt 


-4.53 


N 


SN 2006bt 


2.27 


N 


SN 


2006bz 


-2.44 


N 


SN 


2006cj 


3.43 


N 


SN 2006cm 


-1.15 


? 


SN 


2006cp 


-5.30 


N 


SN 


2006cq 


2.00 


N 


SN 2006cs 


2.28 


F 


SN 


2006CZ 


1.12 


? 


SN 


2006ef 


3.20 


N 


SN 2006gr 


-8.70 


? 


SN 


2006ej 


-3.70 


N 


SN 


2006et 


3.29 


N 


SN 2006gt 


3.08 


N 


SN 


2006ke 


2.36 


N 


SN 


2006kf 


-8.96 


? 


SN 2006kf 


-3.05 


N 


SN 


200611 


-6.30 


? 


SN 


20061e 


-8.69 


N 


SN 200finr 


—2.79 


N 


SN 


2006sr 


— 2.34 


N 


SN 


2006sr 


2.69 


N 


SN 2007A 


2.37 


N 


SN 


2007F 


-9.35 


A 


SN 


2007F 


3.23 


N 


SN 2007N 


0.44 


N 


SN 


2007O 


-0.33 


N 


SN 


2007af 


-1.25 


N 


SN 2007af 


2.84 


N 


SN 


2007al 


3.39 


? 


SN 


2007bc 


0.61 


F 


SN 2007bd 


-5.79 


N 


SN 


2007bm 


-7.79 


A 


SN 


2007bz 


1.65 


N 


SN 2007ca 


-11.14 


N 


SN 


2007ci 


-6.57 


N 


SN 


2007ci 


-1.71 


N 


SN 2007CO 


-4.09 


N 


SN 


2007CO 


0.85 


N 


SN 


2007cq 


-5.82 


A 


SN 2007fb 


1.95 


N 


SN 


2007fr 


-5.83 


N 


SN 


2007fr 


-1.25 


? 


SN 2007gi 


-7.31 


N 


SN 


2007gi 


-0.35 


N 


SN 


2007gk 


-1.72 


N 


SN 2007hj 


-1.23 


N 


SN 


20071e 


-10.31 


N 


SN 


20071c 


-9.40 


N 



Continued on Next Page. . . 
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Table 2 — Continued 



SN Name 


Phase^ 


C II 

Type*' 


SN Name 


Phase'' 


C II 

Type'' 


SN Name 


Phase^ 


C II 

Type'^ 


SN 200781== 


-1.23 


N 


SN 2007on 


-3.01 


F 


SN 2007on 


-3.00 


F 


SN 2007qe 


-6.54 


N 


SN 2008Z 


-2.29 


F 


SN 2008ar 


-8.87 


N 


SN 2008ar 


2.83 


N 


SN 2008bt 


-1.08 


? 


SN 2008S1'' 


-6.36 


A 


SN 2008S1'' 


-4.40 


A 


SN 2008S1'' 


-3.42 


F 


SN 2008S1'' 


0.49 


F 


SN 2008dx 


2.46 


? 


SN 2008ec 


-0.24 


N 


SN 2008ei 


3.29 


N 


SN 2008s5° 


1.26 


N 


SN 2008hs 


-7.94 


F 









''Phases of spectra are in rest-frame days relative to B-band maximum using the heliocentric redshift 
and photometry reference presented in Table 1 of BSNIP I. 



''Classification based on the presence or absence of C II absorption in our spectra. 'A' = clearly 
separated spectral absorption feature attributed to C II is present; 'F' = a flattening or depression 
in the red wing of the Si II A6355 feature is present that is likely due to C II; 'N' = no evidence 
for C II is present; '?' = no deflnitive determination can be made about the presence or absence of 
C II owing to noisy data. 
"Also known as SNF20071021-000. 
''Also known as SNF20080514-002. 
'^Also known as SNF20080909-030. 
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Table 3. Total Number of Objects and Spectra 



C Class. 


# Objects 


# Spectra 


A 


16'' 


19 


F 


20 


29 


N 


95" 


117 


? 


13 


23 


Total 


144 


188 



an additional ~14 per cent show some evidence for C in their 
spectra ('F'), for a total of ~25 per cent 'A' + 'F.' This is 
consistent with previous studies that find ~ 20-33 per cent 



''Classification based on the presence or absence of 
C II absorption in our spectra. 'A' = clearly sepa- 
rated spectral absorption feature attributed to C II 
is present; 'F' = a flattening or depression in the 
red wing of the Si II A6355 feature is present that 
is likely due to C II; 'N' = no evidence for C II 
is present; '?' = no definitive determination can be 
made about the presence or absence of C II due to 
noisy data. 

''Two of these objects were classified as 'N' using 
the BSNIP data, but reclassified as 'A' based on 
data presented in previous work; see text. 
'^Two of these objects were classified as '?' using 
the BSNIP data, but reclassified as 'N' based on 
data presented in previous work; see text. 



There are 34 SNe la in the current sample for which 
we inspect multiple spectra, and in many of these cases the 
C classifications of different spectra of the same object dis- 
agree. However, this is unsurprising since C features tend 
to weaken with time (e.g., |Folatelli et al.||2012"| ). In all cases 
using the BSNIP data, the temporal evolution of the C clas- 
sification is 'A'— ^'F'— >-'N.' Therefore, we classify a SN by 
the C classification of its earliest spectrum. Also, since any 
spectrum classified as '?' does not indicate whether C is 
present, we ignore any '?' spectra when determining the C 
classification of a given SN. 

When comparing the BSNIP data to previous stud- 



ies, 33 SNe la have been classified in earlier work (Parrent 
|et al.|[20TT| [Thomas et"ar][20TT] [Folatelh et"aI1|2012[ ) and 

of these, our classifications agree for 23 objects. For most 
of the objects where the C classification differs, the BSNIP 
spectra are from earlier epochs or have higher S/N and show 
stronger evidence for C than previously thought; for these 
objects we retain the C classification determined from the 
BSNIP data. However, there are two objects we classify as '?' 
that were previously classified as 'N' using higher S/N data 
from earlier epochs, as well as two objects that we classify as 
'N' that were previously classified as 'A' using spectra from 
earlier epochs. In these four cases we adopt the classification 
from the literature in lieu of the one determined from our 
own data. These reclassifications are reflected throughout 
this work. We note here that if we had access to more high- 
quality data or more spectra at younger epochs, we might 
classify even more objects as C positive. Therefore, the in- 
cidence rates of C in SNe la that we calculate herein should 
be considered a lower limit. 

The final C classification for each SN la in our sample 
(including the above reclassifications) can be found in Ta- 
ble [TJ and a summary of the number of objects in each class 
is presented in Table[3] In the BSNIP data, ~11 per cent of 
the SNe la show definite C absorption features ('A'), while 



of SNe la show evidence for C ( 


Parrent et al.|2011 Thomas 


et al.|2011 Folatelh et al.|2012 


1. 



3.2 Measuring the Carbon 

Figure [2] shows the region near C II A6580 and Si II A6355 
of all 19 spectra in our sample that are classified as 'A.' The 
wavelengths corresponding to C II A6580 with expansion ve- 
locities of 10,500-14,000 km s"^ (i.e., the range of velocities 
observed in this work) are highlighted. For each of these 
spectra, we determine the pEW and expansion velocity of 
the C II A6580 feature, and we attempt to measure these pa- 
rameters for the C II A7234 feature as well. The algorithm 
used to measure the C II A6580 absorption is described in 
detail in BSNIP II, but here we give a brief summary of the 
procedure. 

Each spectrum first has its host-galaxy recession ve- 
locity removed and is corrected for Galactic reddening (ac- 
cording to the values presented in Table 1 of BSNIP I), and 
then is smoothed using a Savitzky-Golay smoothing filter 
(Savitzky & Golay 1964). We attempt to define a pseudo- 



continuum for each spectral feature. This is done by deter- 
mining where the local slope changes sign on either side of 
the feature's minimum. Quadratic functions are fit to each 
of these endpoints, and the peaks of the parabolas (assum- 
ing that they are both concave downward) are used as the 
endpoints of the feature; they are then connected with a 
line to define the pseudo-continuum. This defines the pEW 
(e.g., Garavini et al.|2007 l. Once a pseudo-continuum is cal- 
culated, a cubic spline is fit to the smoothed data between 
the endpoints of the spectral feature. The expansion veloc- 
ity is calculated from the wavelength at which the spline fit 
reaches its minimum. Every fit is visually inspected, and the 
fits to all of the 19 'A' spectra are found to be acceptable. 

We attempt to use the same procedure for the C II 
A7234 absorption feature, but due to its relative shallow- 
ness and the difficulty in automatically defining the end- 
points, our fitting routine fails on nearly all of the 'A' spec- 
tra. Therefore, a manual version of the algorithm is used 
(i.e., the endpoints are defined by hand). Even so, the A7234 
absorption cannot be accurately measured in 7 of 19 'A' 
spectra. Also, as a sanity check, this manual version of the 
fitting procedure is used to measure the A6580 feature in all 
of the 'A' spectra and the results are compared to those of 
the more robust, automated fitting method described above. 
All values of pEW and velocity from these two methods are 
consistent within the measured uncertainties, adding cred- 
ibility to the values measured for the A7234 feature. Note 
that throughout the analysis presented here we use only 
pEW and velocity measurements for the A6580 feature as 
determined by our automated fitting routine. The results of 
these measurements, for both C II features inspected, can 
be found in Table ID 
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Figure 2. All 19 of our 'A' spectra showing the region near C II A6580 and Si II A6355. The grey regions cover wavelengths corresponding 
to C II A6580 with expansion velocities of 10,500-14,000 km s~l (i.e., the range of velocities observed in this work with SNe 2002cr 
and 1994D having the smallest and largest velocities, respectively). Each spectrum is labeled with the SN name and rest-frame age. The 
data are all deredshifted and dereddened using the redshift and reddening values presented in Table 1 of BSNIP I, assuming that the 
extinction follows the [Cardelli et al.| ( |1989[ l extinction law modified by |0'Donnell| l |1994[ | . 
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Table 4. Measured Values for C II Lines 



SN Name Phase'' A6580 pEW'' A6580 A7234 pEW'' A7234 



SN 


1994D 


— 12.31 


12.22 


(0.09) 


0.5 


(0.1) 


11.23 


(0 11) 


10.2 


(0 5) 


SN 


1994D 


— 11.31 


13.88 


'0.10) 


2.6 


(0.3) 


11.44 


(DM) 


10.5 


CD ",) 


SN 


1994D 


—9.32 


13.03 


0.09) 


1.4 


(0.2) 


10.11 


(0 10~1 


6.9 


CD 3") 


SN 


1994D 


—7.67 


12.65 


(0.09) 


0.9 


(0.1) 


8.91 


(0.09) 


6.5 


(0 3) 


SN 


1994D 


—5.32 


12.27 


'0.09) 


0.6 


(0.1) 










SN 


1998dm 


-12.48 


13.49 


(0.10) 


1.7 


(0.2) 


11.27 


(0.11) 


6.2 


(0.3) 


SN 


2002cr 


-6.78 


10.72 


(0.16) 


2.1 


(0.3) 


9.60 


(0.10) 


6.0 


(0.3) 


SN 


2002hw 


-6.27 


12.14 


(0.14) 


1.0 


(0.2) 










SN 


2003kf 


-7.50 


12.72 


(0.16) 


0.8 


(0.1) 










SN 


2004ey 


-7.58 


12.94 


(0.09) 


0.2 


(0.1) 










SN 


2005cf 


-10.94 


11.88 


(0.09) 


1.0 


(0.2) 










SN 


2005el 


-6.70 


12.50 


(0.09) 


3.1 


(0.4) 


10.24 


(0.10) 


7.2 


(0.4) 


SN 


2005eu 


-9.06 


12.88 


(0.37) 


0.8 


(0.1) 










SN 


2005iq 


-5.86 


12.02 


(0.17) 


2.7 


(0.5) 


11.31 


(0.11) 


4.1 


(0.2) 


SN 


2007F 


-9.35 


12.73 


(0.16) 


1.1 


(0.2) 


11.01 


(0.11) 


4.3 


(0.2) 


SN 


2007bm 


-7.79 


11.25 


(0.09) 


1.5 


(0.2) 


9.47 


(0.09) 


6.8 


(0.3) 


SN 


2007cq 


-5.82 


11.82 


(0.17) 


1.9 


(0.3) 


9.47 


(0.09) 


4.7 


(0.2) 


SN 


2008sld 


-6.36 


11.82 


(0.16) 


1.9 


(0.3) 


11.18 


(0.11) 


4.0 


(0.2) 


SN 


2008sl<^ 


-4.40 


11.00 


(0.09) 


1.7 


(0.4) 











la uncertainties for each measured value are given in parentheses. 

''Phases of spectra are in rest-frame days using the heliocentric redshift and 

photometry reference presented in Table 1 of BSNIP I. 

''The pEW is in units of A. 

"^The expansion velocity is in units of 1000 km s^^. 
^Also known as SNF20080514-002. 



4 ANALYSIS 

4.1 When is Carbon Detectable? 

As stated above, the oldest 'A' spectrum in the BSNIP sam- 
ple is from 4.4 d before maximum brightness and the old- 
est 'F' spectrum was obtained 3.6 d after maximum. This 
matches well with previous studies, which found that 'A' 
spectra are found at ages less than 3 d before maximum 
(Parrent et al.||2011| iThomas et al.||2011| IFolatelli et al 



2012). Furthermore, all three of the earlier investigations 



only use pre-maximum spectra (for SNe la that are not pos- 
sible super-Chandrasekhar-mass or SN 2002cx-like objects), 
and like the current study they find 'F' spectra near maxi- 
mum brightness. 

The top panel of Figure [3] shows the fraction of spectra 
with C (just A,' and the sum of A' and 'F') as a function of 
time. The horizontal error bars represent the width of each 
bin (i.e., 2 d) and the vertical error bars represent the range 
of fractions if one SN with C is added to or subtracted from 
that bin. The fraction of A' spectra and A'-|-'F' spectra 
both start at 60 per cent at 12 d before maximum bright- 
ness. The fraction of A' spectra decreases monotonically 
with time, while the fraction of A'+'F' spectra generally 
(but not always) decreases with time as well. 

The small rise in the fraction of A'+'F' spectra from 
t = — 12 d to t = — 10 d appears to be insignificant. However, 
the rise of A'-|-'F' spectra for —8 d ^ t ^ —4 d seems in- 
consistent with an actual decrease. It is unclear what might 
cause this small spike in the fraction of spectra with C at 
these epochs. This is consistent with what was seen by |Fo-| 



latelli et al. ( 2012 1 in their Fig. 11, but not exactly the same. 



Their fractions of A' and A'+'F' both monotonically de- 
crease with time, and they detect no 'F' spectra older than 
2 d before maximum. Using the BSNIP data, 10-20 per cent 
of the spectra with —2 d < t < 5 d are classified as 'F.' 

The cumulative fraction of spectra with C (both 'A' and 
'A'+'F') is shown as a function of time in the bottom panel 
of Figure [3] Each point represents the fraction of spectra 
with C from epochs in that bin or younger. The symbols 
and error bars have the same meanings as in the top panel. 
Once again there is a monotonic decrease in the fraction 
of 'A' spectra with time, in addition to mostly decreasing 
fractions of C-positive spectra (i.e., 'A'+'F') with time. By 
t ~ 5 d, which corresponds to the age bin that includes 
our oldest 'F' spectrum, ~12 per cent of the spectra in the 
BSNIP sample show definitive C signatures (i.e., 'A') while 
~29 per cent of them show at least possible evidence for C 
(i.e., 'A'+'F'). 

It seems that if one wants to detect C in an optical 
spectrum of a SN la that follows the Phillips relation, a rel- 
atively high-quality spectrum must be obtained at an epoch 
younger than ~4 d past maximum brightness. However, ob- 
servations at the end of this epoch range will yield only a 
possible C signature (i.e., an 'F' spectrum) and will occur 
<10 per cent of the time. At t ~ —4 d, the chance of de- 
tecting C goes up significantly. For the BSNIP data there is 
about a 50 per cent chance of obtaining an 'A' or 'F' spec- 
trum at this epoch, though the probability of obtaining an 
'A' spectrum is still only ~8 per cent. Finally, it appears 
that obtaining spectra at t < —5 d yields a relatively good 
chance of showing some sign of C (just over 50 per cent) and 
a better than one-third chance of yielding an 'A' spectrum. 
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and SN 1999aa-like objects ("Ia-99a a," |Li et al.|200"T] [StroL 
ger et al.|20"02}|Garavini et al. 12004] ). See BSNIP I for more 
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Figure 3. {Top) The fraction of spectra with C versus time. The 
sum of 'A' and 'F' are asterisks and 'A' alone are squares. The 
horizontal error bars represent the width of each bin (i.e., 2 d) and 
the vertical error bars represent the range of fractions if one SN 
with C is added to or subtracted from that bin. {Bottom) The 
cumulative fraction of spectra with C versus time. This is the 
fraction of spectra with C in that age bin or younger. Symbols 
and error bars have the same meanings as in the top panel. 



4.2 Carbon and Various Classification Schemes 

In Tablejl] the "Carbon Classification" for each object stud- 
ied here is listed, along with its spectral classification based 
on various other classification methods. The "SNID type" of 
each SN is taken from BSNIP I. The SuperNova IDentifica- 



tion code (SNID; Blondin & Tonry 20071, as implemented 
in BSNIP I, was used to determine the spectroscopic sub- 
type of each SN in the BSNIP sample. SNID compares an 
input spectrum to a library of spectral templates in order 
to determine the most likely spectroscopic subtype. Spectro- 
scopically normal objects are objects classified as "la-norm" 
by SNID. 

The spectroscopically peculiar SNID subtypes used here 
include the often underluminous SN 1991bg-like objects 
("Ia-91bg," e.g., Filippenko et al.|fl992b Leibundgut et al. 
|1993[ ), and the often overluminous SN 1991T-like objects 



information regarding our implementation of SNID and the 
various spectroscopic subtype classifications. If an object has 
a SNID type of simply "la," it means that no definitive sub- 
type could be determined. 

According to Table [l] all 'A' objects are la-norm, with 
the exception of one "la." All 'F' objects are also la-norm, 
except one each of Ia-99aa, Ia-91bg, and la. On the other 
hand, all SNID types are well-represented (relative to their 
overall incidence rate) in the 'N' objects. However, due to 
the relative rarity of the spectroscopically peculiar subtypes, 
one would expect only 1-2 of each of the non-Ia-norm sub- 
types in a sample of 19 objects ( Ganeshalingam et al.|2010| 



Li et al. 2011b I. Thus, it seems that the SNe la showing 



evidence for C in their spectra are spectroscopically normal 
objects when examining their entire optical spectrum (as 
SNID does). Note, though, that there could exist very rare 
cases of spectroscopically peculiar SNe la that follow the 
Phillips relation and that show C features. 

The fourth column of Table [l] presents the "Benetti 
type" of each object, which is based on the velocity gra- 
dient of the Si II A6355 feature ( [Benetti et"aI1 |2005|). The 



high velocity gradient (HVG) group has the largest velocity 
gradients while the low velocity gradient (LVG) group has 
the smallest velocity gradients. The third subclass (FAINT) 
has the lowest expansion velocities, yet moderately large ve- 
locity gradients, and consists of subluminous SNe la with 
the narrowest light curves. As mentioned in BSNIP II, from 
where these classifications are taken, the BSNIP data are 
not well-suited to velocity-gradient measurements since the 
average number of spectra per object is ~2 (see BSNIP I). 
However, we are still able to calculate the velocity gradient 
for a subset of our data. 

Of the SNe la with C and a known Benetti type, four are 
LVG (three of which are 'A') and four are HVG (two of which 
are 'A'). We also find that the actual values of the velocity 
gradient itself are similar for objects with and without C. 



Parrent et al. (2011j find that LVG objects have a greater 



chance of showing C as compared to HVG objects and that 
no HVG SNe show a definitive C signature. They point out 
that this could be an observation bias since HVG objects 
tend to have higher Si II velocities (e.g., 'Hachingc r et al?] 
Wang et al.|2009[ ), increasing the amount of blending 



2006 



between Si II A6355 and C II A6580 and thus making it more 
difficult to detect C. In BSNIP II we show that the one-to- 
one association between HVG and high expansion velocities 
is not as clear as has been assumed previously, which could 
explain how we are able to detect C in some HVG objects. 
However, this possible connection between velocity gradient 
and incidence of C should be explored further in the future 
with datasets that are more suited than BSNIP to velocity- 
gradient calculations. 

The "Branch type" referred to in Table [l] uses pEWs of 
Si II A6355 and Si II A5972 measured near maximum light 
to classify SNe la (Branch et al. 20061. The four groups 



("Ia-91T," e.g., lFilippenko et al.|1992a]lPhillips et al.|1992 | 



they define based on these two pEW values are core nor- 
mal (CN), broad line (BL), cool (CL), and shallow silicon 
(SS). However, they point out that SNe seem to have a con- 
tinuous distribution of pEW values; hence, how the exact 
boundaries are defined is not critical. The Branch-type clas- 
sifications used here can be found in BSNIP II. The majority 
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(63 per cent) of CN objects show evidence for C, while only 
16 per cent of BL objects have C in their spectra. This is 
consistent with the idea mentioned above that it is harder 
to distinguish C absorption in SNe la having high expansion 
velocities such as BL objects ( Parrent et al.pOll Folatelli 
et al.|2012| ). 

Furthermore, only 25 per cent and 18 per cent of SS 
and CL objects show evidence for C, respectively, which has 



been noticed in earlier work (Parrent et al. 2011 Folatelli 



et al.|20T2 i. This has been interpreted as evidence that the 



presence or absence of C depends on the effective tempera- 
ture ( Parrent et al.|2011 1, since the effective temperature is 
directly related to the relative strengths of the Si II A6355 
and Si II A5972 features ( [Nugent et al]|1995[ ). 

The prevalence of CN objects with C, as compared to 
other Branch types, is unsurprising given what was found 
above when discussing SNID types. In BSNIP II it was 
shown that SNID types are often equivalent to the more 
extreme objects in each of the non-CN Branch types. There- 
fore, since nearly all SNe la with C are la-norm, it stands 
to reason that most of them should also be CN. 

The last column of Table[T]lists the "Wang type" of each 
object, taken from BSNIP II, which is determined from the 



Si II A6355 velocity near maximum brightness (Wang et al 



2009). SNe la which are classified as la-norm by SNID and 
have high velocities near maximum (>11,800 km s"'^) are 
considered to be high-velocity (HV) objects. la-norm with 
velocities less than this cutoff are classified as normal (N). 

No HV objects are classified as 'A' and only 2 are clas- 
sified as 'F,' while 21 HV objects show no evidence of C. 
Nearly one-third of normal-velocity objects, however, have 
C. This is consistent with the relative lack of BL objects 
that show C, since both BL and HV SNe have large expan- 



sion velocities and this makes C detection difficult (Parrent 
et al.|2011||Folatelh et"al]|2012[ ). 



4.3 Similarities (and Differences) Between 
Objects With and Without Carbon 

As mentioned above, the Phillips relation correlates the peak 



luminosity of a SN la to its light-curve decline rate ( Phillips 



19931. One way to parametrise this decline rate is to cal- 



culate the difference in magnitudes between maximum and 
fifteen days past maximum in the B band, referred to as 
Ami5(_B). The BSNIP sample has 196 objects for which we 
calculate Ami5(_B) and of those, 28 show C and 67 show 
no C. A histogram of Ami5(_B) values can be found in the 
top panel of Figure |4] The average Ami5(_B) for each of the 
three samples show in the figure (all of BSNIP, with C, and 
without C) are consistent with each other. This was hinted 
at by Folatelli et al. (20121, though they admit that they 



have too few objects to make any robust statistical state- 
ment. There also appears to be no difference in Ami5(i3) 
values when the "with C" sample is subdivided into 'A' and 
'F' objects. 

We also use an alternative parameterization of the 
light-curve width: the xi parameter from SALT2 ( |Guy| 



et al. 20071. The value of xi is in the opposite sense as 



Ami5(-B). Thus, underluminous, narrow, fast-evolving light 
curves have large values of Ami5(_B) but small values of xi. 
There are 335 SNe la in BSNIP that have a SALT2 fit, and 
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Figure 4. (Top) A histogram of Ami5(S) for the entire BSNIP 
dataset (grey), with C (red hashed), and without C (blue hashed). 
There is no significant difference between Ami5(B) values for any 
of the three samples. (Bottom) A histogram of xi using the same 
colours as the top panel. Again, there is no significant difference 
between xi values for any of the samples. 



30 (83) of them exhibit (do not exhibit) C. A fiistogram of 
xi values is shown in the bottom panel of Figure |4] 

Like Ami5(_B), the average xi value for each sample is 
consistent with one another. This is at odds with the finding 
of [Thomas et al.| ( |2011[ ) that SNe la with C have lower xi 
values when compared to SNe without C. Their C-positive 
objects are mostly clustered near xi ~ —2 while our SNe la 
with C have a wide range of xi values, with the average 
being xi ~ —0.74 and the peak of the distribution occurring 
at xi ~ 0. However, the overall distributions of xi values 
are different between [Thomas et al.| ( |2011| and the current 
study; there are a significant number of SNe la with xi < 
—2.5 in the BSNIP sample while there are none shown by 
[Thomas et al] ( [20TT| . 

The colours of SNe la that show or do not show C sig- 
natures can also be investigated. One way to quantify the 
colour of a SN la is to measure the difference between its 
B-band magnitude and V-band magnitude at the time of 
-B-band maximum brightness (referred to in this work as 
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Figure 5. (Top) A histogram of the difference between SN B- 
band magnitude and V-band magnitude at the time of S-band 
maximum brightness (i.e., {B — V)^^^) for the entire BSNIP dat- 
set (grey), with C (red hashed), and without C (blue hashed). 
(Bottom) The cumulative distribution function of (B — V)^^^ 
values using the same colours as the top panel. SNe la with C 
tend to have bluer colours than those without C. 



{B - V)^^^). The BSNIP data contain 190 objects for which 
{B — V)^^^ is measured, 28 of which have C and 67 of which 
do not. While there is no significant difference in (B — V)^^^ 
when the "with C" sample is subdivided into 'A' and 'F' 
objects, there is a difference between the "with C" objects 
and the "without C" objects: SNe la with C are bluer. A 
Kolmogorov-Smirnov (KS) test on these two samples implies 
that they very likely come from different parent populations 
(p « 0.07). The top panel of Figure [5] shows a histogram 
of {B — V)^^^ values for the entire BSNIP dataset, objects 
with C, and those without. The bottom panel shows the 
cumulative distribution function of these three samples. 

This trend can also be seen if one parametrises SN 
colour by the SALT2 c parameter ( |Guy et al.|2007[ ). Of the 
335 objects in BSNIP with good SALT2 fits, 30 have C and 
83 do not. And, as with the {B — V)^^^ values, C-positive 
objects appear to have bluer colours than C negative ones (a 
KS test yields p ~ 0.01), while no significant colour differ- 



Figure 6. (Top) A histogram of SALT2 c for the entire BSNIP 
datset (grey), with C (red hashed), and without C (blue hashed). 
(Bottom) The cumulative distribution function of SALT2 c values 
using the same colours as the top panel. Again, SNe la with C 
tend to have bluer colours than those without C. 



ence is found between 'A' objects and 'F' objects. Figure [6] 
presents the histogram (top panel) and the cumulative dis- 
tribution function (bottom panel) of the SALT2 c values for 
each of the three samples. 

Yet another way to quantify the colour of a SN la is to 
calculate synthetic photometric colours from a spectrum. In 
BSNIP I it was shown that the relative spectrophotometry 
of our data, when compared to the actual light curves, is 
accurate to ^ 0.07 mag across the entire spectrum. There- 
fore, synthetic colours derived from BSNIP spectra should 
be photometrically accurate to about this level. In order to 
determine the synthetic colours of our spectra in this work, 
we follow the procedure from BSNIP I. Simply stated, we 



convolve each spectrum with the Bessell ( 1990 1 filter func- 



tions, which have approximate wavelength ranges of 3400- 
4100, 3700-5500, 4800-6900, 5600-8500, and 7100-9100 A 
for U, B, V, R, and I, respectively. We then calculate the 
U - B, B - V, V - R, R - I, and U -V colours. Most of 
the spectra studied herein fully cover the BVRI bands and 
about half cover the U band as well. Using synthetic colours 
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calculated from our spectra, objects with C signatures once 
again have significantly bluer colours at all epochs for all 
colours calculated (except R — I, where SNe la without C 
tend to have bluer colours than those with C). 

The result found here that SNe la with evidence for C 
tend to have bluer colours confirms the work of previous 
groups ( |Thomas et al.||2011| |FolateUi et al.||2012[ ). A rela- 
tionship between colour and light-curve width was shown 



by Folatelli et al. (20121 to be steeper for objects with C as 



compared to objects without C when only SNe la with little 
intrinsic reddening were considered. The BSNIP data show 
no significant evidence of different light-curve width versus 
colour relationships for objects with or without C, whether 
we use all objects or just those that are unreddened. 
Milne et ah] ( |2010[ ), [Thomas et "aL] (|2011[), and 



fc Brown] (|20T2f pr esent Swift/XJYOT ( [Gehrels eFal 



Milne 



2004 



Homing et al. 2005| photometry of a handful of SNe la. In all 



of these works it is shown that the objects that are relatively 
bright in the near-ultraviolet ("NUV blue"; i.e., those with 
the largest NUV excesses) also show strong evidence for C 
absorption. However, they note that SN 2005cf, an object 
which is clearly C positive, has completely normal colours 
in the Sw ift/VVOT data . All o f the NUV-blue objects pre- 
sented by iThomas et al.l (|2011|) and iMilne & Brownl (12012 1 



that are also studied in this work are found to exhibit C fea- 
tures and we find no evidence for C in three SNe la that are 
NUV-red (according to the Swift/UVOT data). However, 
two objects in BSNIP that are NUV-red appear to have C 
signatures in their spectra: SN 2005cf (as mentioned above) 
and SN 2007cq (the Simft/UVOT data indicate that it is 
NUV-red, even though it has quite blue optical colours — 
its {B ~ V)^^^ and SALT2 c are low, 0.004 mag and 0.024, 
respectively) . 

In summary, SNe la which show evidence for C in 
their pre-maximum spectra have bluer colours in the opti- 
cal bands at all pre-maximum epochs and are almost always 
found to be NUV-blue in space-based UV/optical photom- 
etry. On the other hand, all SNe la which are found to be 
NUV-blue in the Swift/UVOT data are C-positive objects. 
Finally, SNe la which show no evidence for C always have 
redder colours in both the optical and NUV. 
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Figure 7. The temporal evolution of the expansion velocity of 
C II A6580 {top) and C II A7234 (bottom). The different shapes 
correspond to different SNe. The two objects which have multiple 
velocity measurements have their velocities connected with a solid 
line. 



4.4 C II Velocities 

As described in Section[3?2j we measure the expansion veloc- 
ity and pEW of the C II A6580 line in all 19 spectra classified 
as 'A.' We also measure the velocity and pEW of the C II 
A7234 line in 12 of these spectra. The temporal evolution of 
the expansion velocities for both features can be found in 
Figure [7] The difi'erent shapes correspond to different SNe, 
and for the two objects which have multiple velocity mea- 
surements (SN 1994D and SN 2008s][^, their velocities are 
connected with a solid line. 

The range of velocities spanned by both features is rel- 
atively small. The C II A6580 line mostly has velocities 
around 12,000-13,000 km s"\ while the C II A7234 line 
is mainly found between 9500 km s~^ and 11,500 km s~^. 
However, this is a somewhat larger range of A6580 velocities 
than what has been seen in previous work ( [Folatelli et al.| 



Also known as SNF20080514-002. 



20121. The largest C II velocity observed in the BSNIP data 



~ 14,000 km s , which may be caused more by an ob- 
servational bias than a real, physical limit. Both |Parrent| 
et al.|(|201l|) and|FolatelU et al.|([2012| discuss the difliculty 



of measuring C II A6580 with v > 15, 000 km s~ due to 
the fact that at these high velocities the feature becomes 
strongly blended with Si II A6355. 

As seen in Figure [7] the typical C II velocities of all 
objects at a given epoch decrease with time (as has been 
seen before; Folatelli et al.|20"l2 |, though there is a dramatic 
increase in velocity between the first and second epochs of 
SN 1994D. This has not been seen in previous work, likely 
due to the fact that our data were obtained at earlier epochs. 

( 2011| ) measure aU C II A6580 
and they see very little 



Furthermore, Thomas et al 
velocities to 



be ~ 12,000 km s" 
change with time. 

The difference in the typical velocities of the two C II 
features implies that the A7234 line is ~2000 km s~^ slower 
than the A6580 line. There has been little attempt previously 
to determine the velocity of C II A7234 due to its relative 
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Figure 8. (Top) The temporal evolution of the Si II A6355 veloc- 
ities for SNe la with C (red circles), without C (blue squares), and 
the Icr region around the average velocity as determined by the 
entire BSNIP sample (grey area). (Bottom) The temporal evo- 
lution of the ratio of the C II A6580 velocity to the Si II A6355 
velocity. The plot symbols are the same as in Figure[7] The dashed 
line is the median ratio (~1.05) and the dotted lines are the me- 
dian ±10 per cent. In both panels, objects with multiple velocity 
measurements are connected with a solid line. 



weakness, but |Thomas et ar] ( |2011[ ) do mention possible de- 
tections of this absorption at velocities somewhat lower than 
those of C II A6580 (consistent with what is found here). 

The Si II A6355 velocities (as measured in BSNIP II) of 
the objects with and without C can also be compared and 
are shown in the top panel of FigurejS] Red circles are SNe la 
with C and blue squares do not have C. The grey shaded 
area represents the Icr region around the average Si II A6355 
velocity from the entire BSNIP sample. Individual objects 
with multiple velocity measurements are connected with a 
solid line. The panel includes all 131 objects in this work 
with a definitive C classification ('A,' 'F,' or 'N'). 

According to the figure, SNe la with C tend to have 
lower Si II A6355 velocities. Nearly all of the C-positive spec- 
tra have Si II velocities that are at or below average, while 
the objects without C span the entire range of Si II veloci- 
ties from below average to well above average. This lack of 



HV objects that show C was mentioned above, and is likely 
in part due to the difficulty in detecting and measuring C II 
A6580 at large expansion velocities. In fact, |Folatelli et al.| 
( [2012[ ) found that for all C-positive objects in their sample, 
the Si II A6355 velocities were <12,500 km s~^. We find six 
spectra of C-positive SNe la that have Si II velocities above 
this value, and three of them are from epochs earlier than the 
earliest ones studied by Folatelli et al. |2012J , when one ex- 
pects even larger expansion velocities for all elements. Thus, 
our findings appear to be consistent with those of previous 
work that the Si II A6355 velocities of C-positive SNe la are 
significantly lower than average. 

The typical Si II A6355 velocities for the C-positive ob- 
jects are 10,000-12,000 km s"'^, very close to the typical C II 
A7234 velocities. However, this is 1000-2000 km s~^ slower 
than the typical C II A6580 velocities, as was also found pre- 
viously ( Folatelli et al.|20"l2 l. The bottom panel of Figure [s] 
shows the ratio of the C II A6580 velocity to the Si II A6355 
velocity for all 19 'A' spectra. The plot symbols are the same 
as in Figure [7[ and again the two objects having multiple 
velocity measurements are connected with a solid line. The 
dashed line is the median ratio (~1.05) and the dotted lines 
are the median ±10 per cent. 

The ratio of these two velocities is remarkably constant, 
especially for t > — 10 d. A similar trend was found by |Par-| 
rent et al. (20111, though their average ratio was slightly 



larger than ours (1.1). For a given object, the ratio may in- 
crease somewhat with time, but with only two objects with 
multiple velocity measurements in our sample it is difficult 
to make any definitive statement about this. However, the 



data presented by Parrent et al. (2011 1 seem to support this 
conclusion as well. Furthermore, we note that the C II ve- 
locities are similar to or larger than the Si II velocities which 
supports the idea of the layered structure of SN la ejecta, 
with unburned C found in layers that are further out (i.e., 
faster expanding) than the layers containing newly synthe- 
sized Si. 

For the spectra with t < — 10 d, the two outliers on the 
low end (the first epoch of SN 1994D and SN 2005cf) both 
have relatively low C II A6580 velocities and higher than 



average Si II A6355 velocities (leading to a small ratio). Par- 
|rent et aT] ( |2011| found no normal SNe la with a ratio much 
less than 1, but their ratio for SN 1994D at t —11 d is 
~1, which matches very well our second epoch of SN 1994D. 
Thus, these abnormally low velocity ratios at early epochs 
may be real and should be investigated further in the fu- 
ture with more early-time spectra. SN 1998dm appears to 
be somewhat of an outlier (<10 per cent above the median 
ratio) at the high end at early times; it has a higher than 
normal C II A6580 velocity with a relatively normal Si II 
A6355 velocity (yielding a larger ratio). 

We also investigated any possible correlations between 
C II velocities and photometric parameters. No significant 
correlations were found between C II velocity and light- 
curve width (parametrised by Ami5(_B) or SALT2 xi) or 
SN colour (parametrised by (B - V)^^^ or SALT2 c). 

4.5 C II pEWs 

The temporal evolution of the pEWs for both features can 
be found in Figure [9] The plot symbols are the same as in 
FigurejT] and the two objects which have multiple pEW mea- 
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Figure 9. The temporal evolution of the pEW of C II A6580 (top) 
and C II A7234 (bottom). The plot symbols are the same as in Fig- 
ure [7] In both panels, objects with multiple pEW measurements 
are connected with a solid line. 



surements are connected with a solid line. The C II A6580 
feature has pEWs which are all <3.5 A, which is consistent 
with what has been seen previously ( jFolatellFet al. 2012j . 
The C II A7234 feature, on the other hand, has not been 
measured before, and we find a range of pEW values of ~4- 

11 A. 

As seen in Section |4.1| the probability of detecting C 
decreases with time, mostly due to a weakening of the C II 
absorption features. Figure [9] shows evidence to support the 
idea that, for the most part, the pEWs of the two C II 
features decrease with time. However, we must point out 
that there are only two objects with multiple pEW mea- 
surements. Interestingly, there appears to be an increase in 
the pEW for one of these objects (SN 1994D) at the earliest 
epochs (—13 d <t < —11 d). While the increase is marginal 
and consistent with no change in pEW for the C II A7234 
feature, the pEW increase is quite significant for the C II 
A6580 feature (this can be seen visually in the top-left panel 
of Figure [2|. 

It was suggested by |Folatelli et al.] ( |2012[ | that one might 
observe such an increase in pEW between 13 and lid before 



maximum brightness, but their data did not extend to suffi- 
ciently early epochs to investigate this further. The expected 
increase in pEW was based on synthetic spectra created us- 
ing a Monte C arlo code ( |Mazzali fc Lucy||1993| |Lucy||1999 



Mazzali|2000 1 as implemented in the analysis of SN 2003du 



(Tanaka et al. 20111. FolatelU et al. (20121 find that the 



synthetic spectra at these early epochs show a strong, red 
emission component of the Si II A6355 feature which tends 
to "fill in" some of the C II A6580 absorption, thus lead- 
ing to a low pEW measurement for C. Furthermore, they 
point out that at these epochs some C is below the pho- 
tosphere (leading to a lower measured pEW), and that the 
large expansion velocities at these times make C II A6580 
become increasingly blended with Si II A6355 (again leading 
to difficulty in measuring pEWs of C II). 



The model used for SN 2003du by [Tanaka etaL] pOTT \ 
required 6.8 x 10"^ Mq of C (mass fraction X (C) = 0.002) 
in the velocity range 10,500 < w < 15,000 km s"\ and the 
pEWs measured from these synthetic spectra are plotted as 



open red squares in Figure 9 of Folatelli et al. (20121. This 
velocity range is consistent with the velocities we find for 
the C II A6580 feature; more impressively, the theoretical 



pEW values shown by Folatelli et al. (20121 are excellent 
matches to the pEWs we measure for SN 1994D. Further- 
more, Folatelli et al. (20121 discuss two other models where 
the amount of C is increased and decreased by a factor of 
four from the SN 2003du value. They state that this range of 
C mass (1.7 x 10^^-2.7 x 10~^ M©) includes objects where 
no C is detected as well as objects that have the largest 
C II A6580 pEWs (at ~1 week before maximum light). The 
pEWs measured from the BSNIP data span a similar range 



of values as the data studied by Folatelli et al. (20121, and 
so we find that this mass range for C also encompasses all 
of our data. 

We plot the temporal evolution of the Si II A6355 pEWs 
(as measured in BSNIP II) of SNe la with and without C 
in Figure [To] As before, red circles are SNe la with C, blue 
squares do not have C, and the grey area is the la region 
around the average Si II A6355 pEW from the entire BSNIP 
sample. Again, individual objects with multiple measure- 
ments are connected with a solid line. Objects which show 
evidence for C tend to have slightly below average Si II 
A6355 pEWs, while objects without C follow the average 
pEW distribution quite well. However, the significance of 
this difference between C-positive objects and C-negative 
objects (or the entire BSNIP sample) is relatively weak. 
This may yet again be the observational bias that as the 
Si II A6355 pEW increases, it becomes more blended with 
the C II A6580 feature (making C detection more difficult). 

As with the C II velocities, no significant correla- 
tions were found between C II pEW and light-curve width 
(parametrised by Ami5(_B) or SALT2 xi) or SN colour 
(parametrised by {B — V)^^^ or SALT2 c). Furthermore, 
no correlations were seen between pEW and synthetic pho- 
tometric colours as derived from the spectra themselves 
(Section 4.3 1. Various spectroscopic luminosity and colour 



indicators (which are defined and discussed at length in 
BSNIP II and BSNIP HI) were also found to be uncorre- 
lated with C II pEW. Moreover, no relationship was found 
between the pEW and velocity of either C II feature. 

The pEWs of other spectral features seen in near- 
maximum spectra of SNe la showed no significant correla- 
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Figure 10. The temporal evolution of the Si II A6355 pEW for 
SNe la with C (red circles), without C (blue squares), and the 
1(T region around the average velocity as determined by the en- 
tire BSNIP sample (grey area). Objects with multiple velocity 
measurements are connected with a solid line. 



tion with C II pEW, except for the so-called Mg II and Fe II 
complexes (both of which have correlation coefficients of ap- 
proximately —0.65 to —0.55 for 11 objects). These blends of 
many IGE spectral lines are defined in BSNIP II and encom- 
pass broad, complex absorption features at 4100-4500 A and 
4500-5200 A, respectively. In BSNIP III it was shown that 
the pEWs of the Mg II and Fe II complexes were relatively 
good proxies for colour, with increased pEW implying a red- 
der colour. The anti-correlation between the pEW of these 
features and the pEW of C II indicates that increased pEW 
of C II implies a bluer colour. 

It has already been shown in this work that C-positive 
SNe la tend to have bluer colours than objects without C 
(Section |4.3[ ). In light of the relationship between the pEW 
of C II and the pEW of Mg II and Fe II, perhaps this re- 
lationship with colour extends further than a simple binary 
splitting of objects with C versus those without C. It seems 
possible that the strength of the C feature is directly related 
to the colour of the SN. Objects where no C is detectable 
have the reddest colours, while objects with some C (i.e., 
low pEWs) have moderate colours, and finally objects with 
the most C (i.e., high pEWs) have the bluest colours. Un- 
fortunately, as mentioned above, the pEW of the C II A6580 
feature is not significantly correlated with any direct mea- 
sure of the SN colour (via light curves or synthetic photom- 
etry from the spectra). The fact that the C II pEW appears 
well-correlated with two pEW-based spectral indicators of 
colour is intriguing, though, and should most definitely be 
investigated further in future studies. 



5 CONCLUSIONS 

In this work we have searched for signatures of unburned C 
from the progenitor WD using a subset of the BSNIP spec- 
troscopic sample. We classify 188 spectra of 144 SNe la with 
ages <3.6 d after maximum brightness as either showing def- 
inite C II absorption ('A'), possibly showing evidence for C 



('F'), definitely not showing C ('N'), or inconclusive ('?'). 
The spectrum-synthesis code SYNOW was used to accurately 
classify all spectra that showed possible evidence for C. The 
primary evidence for C is a distinct absorption line associ- 
ated with C II A6580, though absorption from C II A7234 is 
also sometimes detected. 

We find that ~11 per cent of the SNe studied show 
definite C absorption features, while a total of ~25 per cent 
show at least some evidence for C in their spectra, consistent 
with previous work ( Parrent et al.|20lT Thomas et al.|20"TT 



Folatelli et al. 12012 \ . The detection rate of C decreases with 
time, though C can sometimes be seen at all ages younger 
than ~4 d past maximum brightness. Near 4 d before maxi- 
mum brightness, there is a 50 per cent probability of detect- 
ing C, according to the BSNIP data. If one obtains a spec- 
trum at i < —5 d, then there is a better than 30 per cent 
chance of detecting a distinct absorption feature from C II. 

Nearly all objects that show C are spectroscopically 
normal (as defined by various classification schemes) while 
SNe la without C detections are from all spectroscopic sub- 
types. The velocity gradients of objects with and without 
C have a similar average and range, and C detections and 
velocity gradients do not seem to be related in any way. How- 
ever, we again point out that the BSNIP dataset is not well- 
suited to velocity-gradient measurements. The light curves 
of SNe la with and without C also appear to have the same 
distribution. On the other hand, confirming previous work 



(Thomas et al. 2011 Folatelli et al. 2012), objects with C 



tend to have bluer optical colours than those without, and 
some (but not all) also have strong NUV excesses. This is 
shown with the BSNIP data using a variety of optical colour 
measurements. 

The typical expansion velocity of the C II A6580 feature 
is 12,000-13,000 km s~^, which is somewhat faster than the 
usual velocity measured for the C II A7234 feature (and we 
are the first to carefully study the velocity of this feature). 
The Si II A6355 velocities, measured in BSNIP II, tend to 
be lower than average for C-positive objects, while SNe la 
without C have a wide range of Si II velocities. The ratio of 
the C II A6580 to Si II A6355 velocities is remarkably con- 
stant with time and among different objects, with a median 
value of ~1.05 (consistent with what was found by [Parrent | 
et al |2011| ). 



The pEWs of the C II A6580 and C II A7234 features are 
found mostly to decrease with time, though there is a signif- 
icant increase between ~13 and lid before maximum light. 
This is consistent with the predictions made by [Folatelli] 
et al. (20121 from spectral models based on those presented 



by Tanaka et al. (2011 1. The range of pEWs measured from 



the BSNIP data is consistent with earlier work and impli es a 
range of C mass in SN la ejecta of 2 x 10"^-3 x 10"^ M© (|Fb~ 



[latelli et al.||2012j ). C-positive objects tend to have slightly 
lower than average Si II A6355 pEWs, but at a relatively 
low significance. The pEWs of the Mg II and Fe II com- 
plexes are found to be anti-correlated with the pEW of C II 
A6580, implying that bluer objects should have larger C II 
pEWs. However, this correlation is not seen when comparing 
the pEW of C II to direct measures of SN colour. 

Even though this is the largest set of SNe la for which 
C has ever been searched, there are still only a handful of 
strong C detections and measurements of C absorption fea- 
tures. Other studies using independent datasets (Thomas 
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et al. 2011[ Folatelli et al. 12012 \ and literature searches ( Par- 
rent et al.||2011 1 have also been conducted, and we confirm 



most of their findings at higher significance. Still, many more 
moderate-to- high S/N SN la spectra at early epochs are 
needed to better investigate C and further probe WD pro- 
genitor models and SN la explosion mechanisms. New, large- 



scale transient searches such as Pan-STARRS (Kaiser et al. 



2009 



2002[) and the Paloma r Transient Factory (PTF; Rau et al. 
Law et aL]|2009 1 will be critical to moving this topic 



forward as they find progressively more young SNe of all 
types. One success story already is SN 2011fe (PTFllkly), 
discovered only 11 hr after explosion by PTF in MlOl, the 
Pinwheel Galaxy ( [Nugent et al.||2011[ [Li et al.||2011a[ ). The 
search for and measurement of C in the many spectra of it 
obtained at extremely early epochs will further our quest to 
better understand SNe la. 
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